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Abstract 

Context. We investigate iPTF13bvn, a core-collapse (CC) supernova (SN) in the nearby spiral galaxy NGC 5806. This object was 
discovered by the intermediate Palomar Transient Factory (iPTF) very close to the estimated explosion date and was classified as a 
stripped-envelope CC SN, likely of Type lb. Furthermore, a possible progenitor detection in pre-explosion Hubble Space Telescope 
(HST) images was reported, making this the only SN lb with such an identification. Based on the luminosity and color of the progenitor 
candidate, as well as on early-time spectra and photometry of the SN, it was argued that the progenitor candidate is consistent with a 
single, massive Wolf-Rayet (WR) star. 

Aims. We aim to confirm the progenitor detection, to robustly classify the SN using additional spectroscopy, and to investigate if our 
follow-up photometric and spectroscopic data on iPTF13bvn are consistent with a single-star WR progenitor scenario. 

Methods. We present a large set of observational data, consisting of multi-band light curves ( UBVRI , g'r'i'z') and optical spectra. We 
perform standard spectral line analysis to track the evolution of the SN ejecta. We also construct a bolometric light curve and perform 
hydrodynamical calculations to model this light curve to constrain the synthesized radioactive nickel mass and the total ejecta mass 
of the SN. Late-time photometry is analyzed to constrain the amount of oxygen. Furthermore, image registration of pre- and post- 
explosion HST images is performed. 

Results. Our HST astrometry confirms the location of the progenitor candidate of iPTF13bvn, and follow-up spectra securely classify 
this as a SN lb. We use our hydrodynamical model to fit the observed bolometric light curve, estimating the total ejecta mass to be 
1.9 M 0 and the radioactive nickel mass to be 0.05 M Q . The model fit requires the nickel synthesized in the explosion to be highly 
mixed out in the ejecta. We also find that the late-time nebular r'-band luminosity is not consistent with predictions based on the 
expected oxygen nucleosynthesis in very massive stars. 

Conclusions. We find that our bolometric light curve of iPTF13bvn is not consistent with the previously proposed single massive 
WR-star progenitor scenario. The total ejecta mass and, in particular, the late-time oxygen emission are both significantly lower than 
what would be expected from a single WR progenitor with a main-sequence mass of at least 30 M Q . 

Key words, supernovae: general - supernovae: individual: iPTF13bvn 


1. Introduction 

Type Ibc core-collapse supernovae (CC SNe) have either had 
their envelopes stripped of hydrogen (SNe lb) or stripped of 
both hydrogen and helium in the case of SNe Ic (e.g., Filippenko 
1997). The mass loss could either be due to an extensive wind in 
a single massive star (Groh et al. 2014), or due to binary interac- 
tion (Yoon et al. 2010). 

The discovery of iPTF13bvn was made by the interme- 
diate Palomar Transient Factory (iPTF) (Law et al. 2009) in 
the nearby 1 galaxy NGC 5806 on 2013 June 16.24 (UT), just 
0.57 days past the estimated explosion date (JD 2456459.17; 
Cao et al. 2013). Early-time spectra indicated a likely Type 
lb classification. Furthermore, Cao et al. (2013) used adaptive 

1 d = 22.5 Mpc, p = 31.76 mag (Tully et al. 2009). 


optics (AO) images and Hubble Space Telescope (HST) pre- 
explosion archival images to show that there is a possible pro- 
genitor within 80 milliarcsec (mas) of the estimated location 
of the explosion. The luminosity and colors of this progeni- 
tor candidate are consistent with those of a single Wolf-Rayet 
(WR) star. Cao et al. (2013) further argue that the early-time 
light curves (LCs) and spectra in the optical and near-infrared, 
along with the mass-loss rate estimated from radio observations, 
are all consistent with a WR star as the progenitor. 

WR stars are very massive, with zero-age main-sequence 
masses (Mzams) easily surpassing 30 M e (Maeder 1981; 
Crowther 2007). These stars have very strong winds, result- 
ing in high mass-loss rates (M) sometimes even exceeding 
M ~ 1 0 5 M 0 yr _1 . It is believed that the high mass-loss rate can 
cause the entire hydrogen envelope to be expelled before the star 
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undergoes core collapse (Groh et al. 2013b, 2014). The result of 
this would then be a stripped-envelope Type lb supernova. 

Following the discovery and the possible progenitor de- 
tection, Groh et al. (2013a) used their Geneva stellar evolution 
models (Groh et al. 2013b, 2014) in combination with the radia- 
tive transfer code CMFGEN (Hillier & Miller 1998) to model 
iPTF13bvn. They conclude that the possible progenitor candi- 
date detection and the early-time photometry and spectroscopy 
are compatible with a model where the progenitor of iPTF13bvn 
is a single WR star with a main-sequence mass of 32 M Q . 

In this paper, we expand on the discussion of iPTF13bvn. 
In Sect. 2 we describe our follow-up observations and give de- 
tails of the data reduction. Section 3 provides a confirmation of 
the astrometric identification of the HST progenitor candidate. 
In Sect. 4 we present the filtered LCs and describe the construc- 
tion of the bolometric LC of iPTF13bvn. We use semi-analytic 
arguments based on the model of Arnett (1982) and the method- 
ology developed by Cano (2013) to show that the bolometric 
properties of the SN are not consistent with the progenitor being 
very massive (i.e., Mzams > 30 M Q ). In Sect. 5 we report our 
follow-up spectroscopy, confirm the classification of iPTF13bvn 
as a SN lb and provide velocity measurements of the SN ejecta. 
The latter are used in Sect. 6 together with the hydrodynamic al 
model HYDE (Ergon et al. 2014b, in prep.) to further constrain 
the synthesized nickel and ejecta masses of the explosion and the 
helium-core mass of the progenitor. In Sect. 5 we also use late- 
time photometry (> 200 days past the explosion) to constrain 
the amount of oxygen in the ejecta by comparing our data to the 
detailed nebular modeling by Jerkstrand et al. (2014b). These re- 
sults are also inconsistent with a very massive progenitor. 

2. Observations and data reduction 

The discovery of iPTF13bvn was made with the Palomar Oschin 
Schmidt 48-inch (P48) telescope. Cao et al. (2013) also re- 
ported follow-up photometry obtained with the P48, the robotic 
Palomar 60-inch telescope (P60; Cenko et al. 2006), and the 
Las Cumbres Observatory Global Telescope network (LCOGT; 
Brown et al. 2013) up to 20 days past the discovery. Here we 
present photometry from the same telescopes up to 90 days past 
discovery, as well as additional data from the Nordic Optical 
Telescope (NOT) at La Palma obtained around 240 days past 
discovery. 

To produce LCs from the P48, P60, and NOT data, we have 
used an image-subtraction procedure with templates consisting 
of images of NGC 5806 obtained approximately one year before 
the discovery of iPTF13bvn. The point-spread function (PSF) 
is determined and matched prior to image subtraction, and is 
subsequently used for PSF fitting photometry on the subtracted 
frames. For the images obtained by LCOGT we lack templates, 
and thus we estimate the galaxy contribution by fitting a low- 
order surface, and then performing PSF fitting photometry. We 
calibrate our P48 and Sloan filter data against a minimum of 10 
SDSS (Ahn et al. 2014) stars in the field. The Johnson-Cousins 
UBVRI filter data were calibrated against Landolt standard stars 
(Landolt 1992) observed during photometric nights. 

Optical and near-infrared spectra of iPTF13bvn were ob- 
tained with many different telescopes and instruments, start- 
ing within 24 hours after the discovery. Spectra until 16 days 
past the discovery were published by Cao et al. (2013). In 
this paper, we present 6 additional later-time optical spectra 
taken 18-86 days past the SN discovery. These were obtained 
with the Palomar 200-inch telescope (P200) using the Double 
Spectrograph (DBSP; Oke & Gunn 1982), with the NOT using 



Figure 1: Progenitor identification of iPTF13bvn based on reg- 
istering an HST WFC3 image of the SN to a stacked archival 
pre-explosion HST ACS image. The left panels show an image 
taken when iPTF13bvn was clearly visible, and the panels to the 
right show the pre-explosion image. The field of view in the up- 
per panels is 16" x 16". The location of the SN is marked by 
the grey dashed square in the upper panels, and the lower panels 
show this region in detail. Common point sources used for the 
registration are marked with boxed black circles. The intensity 
contours of the bright SN in the lower-left panel are traced with 
black dashed lines. The centroid of the SN is indicated with red 
circles in each panel; the smaller and larger red circles respec- 
tively represent lcr and 5cr uncertainties in the registration. The 
lower-right panel shows the progenitor identification; the cen- 
troid of the SN is directly coincident with one object. 


the Andalucia Faint Object Spectrograph (ALFOSC), and with 
the two Keck 10 m telescopes using the Low Resolution Imaging 
Spectrograph (LRIS; Oke et al. 1995) on Keck 1 and the Deep 
Extragalactic Imaging Multi-Object Spectrograph (DEIMOS; 
Faber et al. 2003) on Keck 2. 

All spectra were reduced using standard pipelines and pro- 
cedures for each telescope and instrument. All spectral data 
and corresponding information is available via WISeREP 2 
(Yaron & Gal- Yam 2012). More data and details will be pub- 
lished by Fremling et al. (2014b, in prep.), where we will also 
investigate the explosion site of iPTF13bvn in more detail. 

3. Progenitor identification 

A progenitor candidate for iPTF13bvn was identified in pre- 
explosion HST images by Cao et al. (2013) 3 . This was done 
by registering AO images taken with OSIRIS and the laser- 
guide-star-AO system on the Keck 1 telescope. Cao et al. (2013) 
present a single source offset by about 80 + 40 mas from the SN 
position, barely within the estimated 2<x image-registration un- 
certainty. We have used recent archival HST (WFC3) images 4 of 

2 http://www.weizmann.ac.il/astrophysics/wiserep/ 

3 First reported in ATels 5140 and 5152 (Arcavi et al. 2013a, b). 

4 Obtained on 2013 Sep. 2.37 (UT), GO-12888, PI S. Van Dyk. 
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Figure 2: Multi-band light curves of iPTF13bvn. Obtained with 
the P48 and P60 telescopes, the LCOGT telescopes, and the 
NOT. 


iPTF13bvn to corroborate this identification. Figure 1 shows the 
result of registering a WFC3 image (filter F555W) to a stacked 
(filters F814W, F555W, and F435W) archival HST (ACS/WFC) 
image from 2005 (Smartt et al. 2009) using 25 common point 
sources. 

Our registration procedure is based on standard MATLAB 
functions and is set up as follows. The centroids of the selected 
common point sources in the images are derived from fits of two- 
dimensional Gaussians that also allow for rotation; the geomet- 
ric transformation is subsequently determined as a second-order 
polynomial transformation with 12 parameters, and the polyno- 
mial transformation is finally applied using bicubic interpola- 
tion. Using this procedure we achieve a standard deviation of 
0.18 HST (ACS/WFC) pixels 5 in the offset between the post- 
explosion centroids and the pre-explosion centroids after image 
registration. If we combine this standard deviation with the un- 
certainty in the centroid of the SN in the post-explosion image, 
we find a total ( 1 <x) uncertainty of 9 mas in our determination of 
the location of iPTF13bvn, corresponding to a projected distance 
of 1 pc. We find that the center of the SN explosion is almost 
exactly coincident (+1.7 mas in a, +3.6 mas in 6) with the cen- 
troid of the previously proposed progenitor candidate. Cao et al. 
(2013) reported the absolute (Vega) magnitudes for this source: 
M b = -5.52 ± 0.39, M v = -5.55 + 0.39, and M, = -5.77 + 0.41 
mag. These are consistent with a single WR star (Eldridge et al. 
2013; Groh et al. 2013a). However, Cao et al. (2013) also cau- 
tion that other scenarios (e.g., a binary system, or a small cluster) 
can result in similar colors and absolute magnitudes. 


4. Light curves 

The P48 r'-band pre-explosion limits and the optical LCs of 
iPTF13bvn are displayed in Fig. 2. We have constructed a quasi- 
bolometric LC 6 by integrating the flux in the BVRI bands, and we 
have determined the black-body (BB) parameters of iPTF13bvn 
by fitting BBs to the spectral energy distribution (SED) derived 
from the full set of UBVRI and g'r'i'z' LCs from the P60 and the 

5 The HST images used here have a pixel scale of 50 mas. 

6 We construct a full bolometric LC and a hydrodynamical model of 
the LC in Sect. 6. 


LCOGT telescopes. For the quasi-bolometric LC calculation, we 
have interpolated the flux in the BVRI bands to the dates of the V- 
band measurements. For these calculations we adopt a distance 
modulus i-i = 31.76 mag (Tully et al. 2009), a Milky Way (MW) 
color excess E(B - V)mw = 0.0278 mag (Schlafly & Finkbeiner 
2011), a host-galaxy color excess E(B - V)h os t = 0.0437 mag 
(Cao et al. 2013), and the Cardelli et al. (1989) extinction law 
with Ry = 3.1. For the BB fits we have also used weights pro- 
portional to the uncertainties in the photometry, and we have 
interpolated or linearly extrapolated the flux in the UBVRI and 
g'r'i'z 1 bands to the dates of both the V- and r'-band measure- 
ments. 

The quasi-bolometric LC and the BB parameters of 
iPTF13bvn are shown in Fig. 3 and Fig. 4, respectively. In 
Fig. 3 we also display the quasi-bolometric BVRI LC of the 
Type lib SN 201 ldh (Ergon et al. 2014a), as well as a scaled and 
stretched quasi-bolometric BVRI LC of SN 1998bw (SN Ic-BL; 
Clocchiatti et al. 201 1). We find that the peak BVRI luminosity 
of iPTF13bvn is somewhat lower than that of SN 201 ldh, in- 
dicating a slightly lower 56 Ni mass 7 . Furthermore, the width of 
the LC peak appears to be narrower for iPTF13bvn than for both 
SNe 1998bw and 201 ldh. The narrow peak already disfavors a 
massive single-star (WR) progenitor scenario, since such a mas- 
sive progenitor (i.e., Mzams ^ 30 M 0 ) should have a large ejecta 
mass resulting in a slower evolution of the bolometric LC with a 
broader peak (Arnett 1982). 

We find that the BVRI LC of SN 1998bw, multiplied by the 
scaling factor k = 0.125 and stretched by the factor ,v = 0.82 
using 

iPTF13bvn LBVRI (f) = k x SN1998bw LBVRI (f x s) (1) 

is very well matched with our BVRI LC of iPTF13bvn. These 
parameter values were found by varying s and k and minimiz- 
ing the residual between the two SN LCs up to 40 days past 
the estimated explosion dates. Interpolating the s and k fac- 
tors found for the sample of SNe Ibc in Cano (2013) provides 8 * 
M56 Ni = 0.05 + 0.02 M 0 and M e j = 1.8 ± 0.3 M 0 . These in- 
terpolations only include SNe with photospheric velocities in 
the range 7000-11000 km s -1 , to roughly match the velocities 
derived from the early-time spectra of iPTF13bvn in Cao et al. 
(2013), and also in later spectra (Sect. 5). A comparison of 
our derived masses for iPTF13bvn to the median nickel mass 
(M.% Ni = 0.15-0.18 M 0 ) and ejecta mass (M e j = 3.9±1.6M 0 )of 
the full SN Ibc sample in Cano (2013) indicates that iPTF13bvn 
has a lower 56 Ni mass and also likely a lower ejecta mass than 
the average SN Ibc. 

From our black-body fits to the SED of iPTF13bvn we de- 
rive the best-fitting BB radius (Rbb) and temperature (7’bb) as a 
function of time. We show the evolution of Rbb, which can be 
interpreted as a rough approximation of the photospheric radius, 
in the bottom panel of Fig. 4. The best-fitting 7’bb is shown in the 
top panel of Fig. 4. Both Rbb and 7’bb are similar to those of SN 
201 ldh (Ergon et al. 2014a), with 7’bb peaking at ~ 8500 K and 
Rbb peaking at ~ 1.4 x 10 15 cm for both objects. Compared to 
the larger sample of SNe Ibc presented by Taddia et al. (2014), 
both 7’bb and Rbb appear to be consistent with those of normal 
SNe Ibc. If Rbb is linearly extrapolated to R B b = 0, we find an 
explosion date (to) of 2013 June 15.75 ± 0.3 (UT), consistent 
with the to estimate of June 15.67 by Cao et al. (2013) using a 
power-law fit to the early-time r'-band LC. 

7 56 Ni mass estimates for SN 201 ldh are 0.075 ± 0.025 M Q by 
Ergon et al. (2014a) and ~ 0.07 M 0 by Shivvers et al. (2013). 

8 Cano (2013) assume a constant opacity, k = 0.07 cm 2 g -1 , in the 

calculations of the physical parameters of the SNe in their sample. 
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Figure 3: Quasi-bolometric BVRI LC of iPTF13bvn (blue cir- 
cles), compared to the quasi-bolometric LCs of SN 201 ldh 
(gray circles) and SN 1998bw (thick dashed line). The LC of 
SN 1998bw has been scaled by A: = 0.125 and stretched by 
s = 0.82. The figure also shows the estimated bolometric LC 
of iPTF13bvn (red circles) compared with our hydrodynamical 
model of the bolometric LC (thin dashed line). 

5. Spectra 

Cao et al. (2013) suggested a Type lb classification of 
iPTF13bvn, partly based on the tentative emergence of Hei 
features (5876, 6678, 7065 A) in optical spectra obtained 10- 
15 days past explosion. Our later-time optical spectra clearly 
demonstrate the continued emergence of these features (Fig. 5), 
as well as of the He i 5016 A line. At ~ 35 days past explosion, 
we find a very good match between the spectra of iPTF13bvn 
and those of the Type lb SNe 2008D (Modjaz et al. 2009) and 
2007 Y (Stritzinger et al. 2009) (bottom panel of Fig. 5). Thus, it 
is clear that the Type lb spectral classification holds. 

SYNOW fits were performed by Cao et al. (2013) on the 
early-time spectra of iPTF13bvn, suggesting a photospheric ve- 
locity (v p h) of 10000 km s _1 at +3 d and 8000 km s -1 at 
+ 12 d. To deduce the photospheric velocity, the absorption min- 
imum of the Fen 5169 A line is likely a more robust tracer 
(Dessart & Hillier 2005). We detect the absorption minimum 
from this line in spectra later than +10 d post-explosion, and 
show the deduced photospheric velocities in Fig. 6. We find 
Vph ~ 10000 km s- 1 at +10 d, in agreement with the 8000- 
10 000 km s' 1 estimate by Cao et al. (2013). Figure 6 also shows 
the photospheric velocity deduced from BB fits to the photome- 
try (Sect. 4), as well as line velocities from the absorption min- 
imum of the He i 5876, 6678, 7065 A lines. Typical uncertain- 
ties in the velocity estimates are 15% before 15 d past the ex- 
plosion and 5-10% at later times. Up until +10 d the differ- 
ent velocity estimates appear to be in good agreement. At later 
stages, especially the Hei 5876 A velocity starts to deviate, and 
ceases to be a useful tracer. Similar behavior was observed for 
SN 201 ldh (Ergon et al. 2014a). In conclusion, we confirm that 
a photospheric velocity of 8000-10000 km s“ during the LC 
peak of iPTF13bvn is a reasonable estimate. For our hydrody- 
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Figure 4: Black-body temperature (top panel) and black-body ra- 
dius (bottom panel) of iPTF13bvn derived from BB fits to the 
photometry. 

namical model of the SN in Sect. 6 we use the Hei 5876 A line 
as a tracer up until +10 d and thereafter the Fe n 5169 A line up 
until +40 d. 

5. 1 . Constraining the oxygen mass 

Our optical spectra at epochs later than +50 d show a clear emer- 
gence of O i features (5577, 6300, 6364, and 7774 A). For the 
massive single-star progenitor scenario proposed by Groh et al. 
(2013a), a large O mass is expected, and our observations allow 
us to test this. 

The concept here is the realization that more massive stars 
produce larger amounts of metals, and in particular the oxygen 
nucleosynthesis is a strong and monotonic function of Mzams 
(Woosley & Heger 2007). This was exploited by Jerkstrand et al. 
(2014a), who find a strong dependence between the ratio (O r ) 
of the [O i] 6300, 6364 A line luminosity to the total 56 Co de- 
cay power, and Mzams of SNe IIP (see, e.g., the top panel of 
Fig. 4 in lerkstrand et al. 2014a). More recently, Jerkstrand et al. 
(2014b) have also modeled nebular stripped-envelope CC SNe, 
and a similar strong dependence between O r and Mzams was 
found. Our spectrum obtained at +86 d is likely not yet in the 
optically thin nebular phase. However, the [Oi] 6300, 6364 A 
line luminosity we measure at this phase is consistent with that 
of SN 201 ldh, which likely was the result of the explosion of 
a 12-13 M e Mzams star with an oxygen mass of 0.3-0. 5 M 0 at 
the time of the explosion (Jerkstrand et al. 2014b) 9 . 

Lacking spectra of iPTF13bvn in the nebular phase, we can 
instead use the r'-band photometry from the NOT obtained at 
+227 d, +240 d, and +262 d. Under the (conservative) assump- 
tion that all of the r'-band flux is due to the [Oi] lines and with 
M.« Ni = 0.05 M 0 , we find O r ~ 0.018 at +227 d, O r « 0.017 
at +240 d and O r ~ 0.015 at +262 d. For SN 201 ldh the cor- 
responding values are in the range O r ~ 0.012-0.08, as derived 
from spectral observations between +200 d to +300 d. For stars 

9 The nebular spectrum of SN 201 ldh has also been modeled by 
Shivvers et al. (2013) who find Mza ms = 13-15 M 0 with an oxygen 
mass of approximately 0.3 M 0 at the time of explosion. 


4 


C. Fremling et al.: The rise and fall of the Type lb supernova iPTF13bvn 



4000 4500 5000 5500 6000 6500 7000 7500 8000 8500 

Rest wavelength (A) 


Figure 5: Spectral evolution of iPTF13bvn in the wavelength range 4000-8000 A (top panel). Comparison of the visible spectra of 
iPTF13bvn to those of other SNe lb at approximately 35 d past the explosion (bottom panel). Thick dashed lines mark the central 
wavelength of the marked emission lines at rest. Thin dashed lines mark the absorption minima associated with the emission lines 
derived from Gaussian fits to the absorption features. The spectra up to 16 d past the explosion have been selected from the early 
spectral sequence presented by Cao et al. (2013). 
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Figure 6: Line velocities of iPTF13bvn derived from absorption 
minima in the spectra (colored circles), the photospheric veloc- 
ity derived from BB fits to the photometry (thick dashed line) 
and the photospheric velocity of iPTF13bvn modeled by our hy- 
drodynamical code (thin dashed line). 


with Mzams significantly higher than that of SN 201 ldh, much 
higher values for O r are predicted; for example, O r ~ 0.035 for 
Mzams = 17 M 0 (Jerkstrand et al. 2014b). Thus, we find that 
the late-time oxygen line luminosity of iPTF13bvn is not con- 
sistent with that expected from the Mzams = 32 M Q progenitor 
suggested by Groh et al. (2013a). 


6. Light curve modeling and progenitor constraints 

The semi-analytic framework of Piro & Nakar (2013), in com- 
bination with early-time observations, was used by Cao et al. 
(2013) to constrain the radius of the progenitor of iPTF13bvn to 
a few solar radii. Furthermore, the HST progenitor candidate de- 
tection and the deduced mass-loss rate by Cao et al. (2013) led to 
the conclusion by both Cao et al. (2013) and Groh et al. (2013a) 
that the progenitor of iPTF13bvn was consistent with a compact 
WR star. However, in the previous sections we have used semi- 
analytic arguments regarding the bolometric LC, as well as neb- 
ular oxygen emission considerations, to show that our post-peak 
data are inconsistent with a massive (> 30 M e ) single-star pro- 
genitor. In this section, we attempt to more robustly constrain 
some of the properties of the progenitor of iPTF13bvn using a 
hydrodynamical model for the bolometric LC. 

We use our observed multi-band UBVRI and g'r'i'z' LCs of 
iPTF13bvn, the explosion date estimate by Cao et al. (2013) as 
an initial guess for the explosion time, and our measured photo- 
spheric velocities together with a grid of SN models constructed 
with the hydrodynamical code HYDE, previously used to model 
the first 100 days of the bolometric LCs of SNe 1993J, 2008ax, 
and 201 ldh (Ergon et al. 2014b, in prep.). The code is one- 
dimensional, based on flux-limited diffusion, and follows the 


framework described by Falk & Arnett (1977) 10 . The model grid 
is based on bare He cores evolved until the verge of core collapse 
using MESA 11 (Paxton et al. 2010). For the infrared and ultra- 
violet regions we use bolometric corrections derived from the 
bolometric LC of SN 201 ldh. When we perform the fits to our 
model grid, we place equal weights on the diffusion phase of the 
LC (+1 d to +40 d), the early tail of the LC (+40 d to +100 d), 
and the photospheric velocity measurements (+2 d to +40 d). 

The best-fitting model, shown in Fig. 3 along with the esti- 
mated bolometric LC, constrains the total energy of the explo- 
sion to be £ = 0.858)'4 qX 10 51 erg, the synthesized radioactive 
nickel mass to be M56 Ni = 0.049^^2 Mo- and the total ejecta 
mass to be M e j = 1.94^° ^ M e under the assumption that a 
1 .5 M 0 remnant remains at the center of the SN explosion. The 
total mass of the He core before the explosion in the model is 
thus Mn e = 3.44+Qgg M 0 . The best-fitting explosion time ( tb ) is 
2013 June 15.55 (UT). The ejecta mass we find is much lower 
than what is suggested by the Groh et al. (2013a) massive single 
star progenitor scenario, which requires 12 M e j » 8-9 M 0 assum- 
ing a 2-3 M g remnant. We also find that to fit the early rise of 
the LC up to 10 d past the explosion, our model requires the ra- 
dioactive nickel synthesized in the explosion to be highly mixed 
out in the ejecta. If the nickel is concentrated toward the center, 
the diffusion time becomes too large, and the rise of the model 
LC happens too late to fit the fast rise of the observed early LC. 
The reported errors are propagated from the observed quantities; 
we assume an uncertainty in the distance modulus to NGC 5806 
of +0.3 mag, an uncertainty in the total extinction for iPTF13bvn 
of ±0.045 mag, and uncertainties in the photospheric velocities 
of +15%. 

We note that our model is generally able to achieve an ex- 
cellent fit to the observed data, including the peak and the tail of 
the LC as well as the photospheric velocities 13 . Furthermore, if 
we let the final explosion time vary between the time of the first 
non-detection until the first detection in the /-band, we find that, 
while the quality of the fit is decreased (especially for dates ear- 
lier than tb - 0.25 d), the high nickel-mixing requirement always 
remains and the derived explosion parameters are not signifi- 
cantly changed. In this sense, our M56 Ni and M e j estimates are 
quite robust, and we note that they are also very similar to the 
semi-analytic nickel and ejecta masses derived in Sect. 4. 

We also note that the radius constraint of a few solar radii 
on the progenitor of iPTF13bvn by Cao et al. (2013) was partly 
based on assumed typical explosion parameters for a SN lb. 
Our model does not directly constrain the radius; it is constant 
across the models in our grid. However, the explosion parame- 
ters we derive permit a larger radius of the progenitor within the 
Piro & Nakar (2013) semi-analytic framework. 


7. Summary and conclusions 

We have confirmed the classification of iPTF13bvn as a Type lb 
SN, and also that the pre-explosion HST progenitor identification 
suggested by Cao et al. (2013) is very likely correct (Sect. 3). 


10 After the submission of this paper, another paper on iPTF13bvn 
(Bersten et al. 2014) has also appeared. Similar hydrodynamical mod- 
eling is performed, showing consistent results. 

11 http://mesa.sourceforge.net 

12 The total He-core mass of the progenitor proposed by Groh et al. 
(2013a) immediately prior to the explosion is 10.9 M 0 . 

13 We do not fit for the photospheric velocities at times later than 
+40 d. Our model of the photosphere is not applicable at these times. 
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While the color and luminosity of the progenitor in the HST im- 
ages are consistent with those of a WR star, we have shown that 
the later-time photometry (Sect. 4) of iPTF13bvn is not consis- 
tent with a single, very massive WR progenitor scenario. Our hy- 
drodynamical model can fit the observed bolometric light curve 
to constrain the total ejecta mass to be ~ 1.9 M 0 . The corre- 
sponding mass of the He core of the progenitor is ~ 3.4 M 0 
immediately prior to the explosion. The synthesized nickel mass 
is constrained to be ~ 0.05 M Q . The model also requires the syn- 
thesized nickel to be highly mixed out in the ejecta, consistent 
with the high average mixing found for the sample of SNe Ibc 
by Taddia et al. (2014). The total ejecta mass and our limits on 
the late-time oxygen emission are both inconsistent with what 
would be expected from a single massive WR progenitor with 
Mzams ~ 30 M 0 as suggested by Groh et al. (2013a). We note 
that while it could be argued that a significant part of the He 
core could fall back onto the central compact object, such a sce- 
nario that still produces an apparently nickel-powered LC, nor- 
mal observed photospheric velocities, and the measured oxygen 
emission appears unlikely. We also note that while the deduced 
high mass-loss rate was interpreted by Groh et al. (2013a) as fa- 
voring a massive WR star, the mass-loss estimate depends on 
an assumed high stellar wind velocity. A lower wind velocity in 
combination with the radio measurements of Cao et al. (2013) 
can also be consistent with a less massive star. 

To date, the SN Ibc showing the largest deviation in LC 
shape compared to the typical SN Ibc LC (e.g., Drout et al. 2011; 
Taddia et al. 2014) is the bolometric LC of SN Ic 201 lbm pre- 
sented by Valenti et al. (2012). SN 201 lbm has a very high peak 
luminosity and (most importantly) a very slow evolution of the 
bolometric LC and a broad peak, which is difficult to explain 
without assuming a very high (> 30 M e ) main-sequence mass 
for the progenitor. This is the kind of SN LC we would expect 
from a scenario including the very massive progenitor star sug- 
gested by Groh et al. (2013a). It is clear that the bolometric LC 
of iPTF13bvn is very different from that of SN 201 lbm. 

If the progenitor is less massive, a binary system for the pro- 
genitor is perhaps the most natural conclusion. Using the pro- 
genitor constraints from our hydrodynamical model, the evolu- 
tionary models for binary systems by Yoon et al. (2010) predict 
values for the radius, mass-loss rate, and hydrogen content in the 
progenitor that are all consistent with what can be derived from 
the early-time observations by Cao et al. (2013). In contrast, cur- 
rent single-star evolutionary synthesis models (e.g., Groh et al. 
2013b) have a harder time producing a progenitor with a suf- 
ficiently low mass to match the observed properties; stars with 
low enough ZAMS masses to explain the low ejecta mass, as 
well as the observed nebular oxygen emission constraints, are 
not predicted to produce SNe Ibc in this context. 

The ultimate test to assess the nature of the progenitor will 
be to reobserve the location of iPTF13bvn with HST after the 
SN has faded. A single massive WR progenitor scenario pre- 
dicts that the progenitor candidate would completely disappear, 
while a less massive progenitor in a binary system would lead to 
a smaller decrease in the luminosity at the location of the pro- 
genitor. 

Finally, we emphasize that the use of the new 
Jerkstrand et al. (2014b) nebular modeling in conjunction 
with photometry during the nebular phase could be an efficient 
way to constrain the ZAMS masses for a larger sample of 
stripped-envelope CC SNe. Future studies of nebular spectra of 
iPTF13bvn are planned, and we are also performing a detailed 
investigation of the host environment of the SN (Fremling et al. 
2014b, in prep.). 
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